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Abstract 

We study analytically the effect of a correlated random potential on the persistent current in a 
one-dimensional ring threaded by a magnetic flux (p, using an Anderson tight-binding model. In 
our model, the system of A^ = 2M atomic sites of the ring is assumed to be partitioned into M pairs 
of nearest-neighbour sites (dimers). While the individual atomic site energies are assumed to be 
identically distributed gaussian variables with autocorrelation parameter Eq, the dimer site energies 
are chosen to be correlated with a gaussian strength a^ < Eq. For this system we obtain the exact 
flux-dependent energy levels to second order in the random site energies, using an earlier exact 
transfer matrix perturbation theory. These results are used to study the mean persistent current 
generated by Ag < A spinless electrons occupying the Ae lowest levels of the flux-dependent energy 
band at zero temperature. Detailed analyzes are carried out in the case of low filling of the energy 
band (1 < Ae < A) and for a half-filled band (Ag = A/2), for magnetic fluxes -1/2 < (p/(f>o < 1/2. 
In the half-filled band case, the uncorrelated part of the disorder reduces the persistent current 
while the correlated part enhances it, in such a way that for a^ < Eq/2 the current decreases with 
the disorder, while for a^ > Eq/2 it increases with it. Also, while showing a specific dependence 
on the flux, the disorder effect has the same dependence on the parity of Ae as the pure system 
free electron current. In contrast, at low filling of the energy-band, the disorder induced effect 
in the persistent current depends critically on the parity: due to a peculiar dependence on the 
flux, it yields a reduction of the current for odd Ag and an enhancement of it for even Ag. The 
observability of the effects of weak correlated disorder on persistent current in the half-filled band 
case is restricted to ring sizes in the nanoscale range, for which no measurements presently exist. 

PACS numbers: 72.15.Rn,73.23-b,73.23.Ra,73.63.Nm 
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I. INTRODUCTION 

The study of persistent current in small metallic rings threaded by a magnetic flux has a 
long history, starting with the basic papers of London[l] and Hund[2]. It has been boosted in 
recent years with the appearance of the seminal paper of Biittiker, Imry and Landauer[3] for 
one-dimensional rings. Important aspects related to the effects of disorder and of electron- 
electron interaction on the persistent current have been reviewed and compared with exper- 
imental observations[4, 5, 6] in well-known monographs[7, 8]and in review articles[9]. An 
important unsolved problem remains however to reconcile theoretical and experimental re- 
sults for the persistent current, particularly in the case where the current averaged over many 
experimental realizations of isolated disordered rings of a given metal has been measured[4]. 
In this case the theoretical results for the ensemble-averaged current are between one or 
two- orders of magnitude lower [7] than the experimental values [4] depending on whether 
electron-electron interaction effects are included or not. However, the electron interaction 
effects are not easy to calculate and no definitive answer as to their precise role seems to 
exist since there are even models in which Coulomb interactions act to reduce the persistent 
current rather than to enhancing it [10]. 

In the present paper we examine another effect which may substantially enhance the 
persistent current in a disordered ring subjected to a magnetic flux. This effect exists when 
the disordered potential is not perfectly random as a result of short-range correlations. 
More precisely, we consider an Anderson tight-binding model of a ring with an even number 
N = 2M of one level atomic sites (of spacing a = 1) whose energies e„ fluctuate randomly 
about a fixed free atom level chosen as zero of energy. As usual the A^ site energies are 
assumed to be identically distributed independent random variables. In order to account 
for short-range correlations we extend the model as follows: we divide the ring up into 
pairs of nearest-neighbour atomic sites (dimers), namely the pairs (1,2), (3,4) ... (2M-1, 
2M), where the sites within individual pairs are assumed to be gaussian correlated, whereas 
sites belonging to different pairs are uncorrelated. Including the usual autocorrelation, the 
Gaussian correlation of site energies is then described by the averages 

(^ra) = 0, \Em£n) = ^^o^min + Oi {Sm,2p-l^n,2p + ^m.,2p^n,2p-l) ■ (1) 

The intersite coupling term in (1) describes short-range structural (dimer) correlations of 



an otherwise uncorrelated (white-noise) random potential. In the following we refer to the 
Gaussian model (1) as an Anderson random dimer model (ARDM) which was first discussed 
in the context of localization on a linear chain[ll], in the case o? = Eq. It differs from 
the random dimer model (RDM) introduced earlier by Phillips and co-workers [12, 13] for 
demonstrating the existence of delocalized states in a one-dimensional system with correlated 
disorder. This model considers chains of A^ lattice sites composed of clusters of host atoms 
a of site energy Sa separated by clusters of defect atoms of energy e^ occurring randomly but 
in pairs of nearest neighbours [12, 13]. Apart from its binary character this system differs 
from our ARDM in that clusters of a-atoms are not restricted to even numbers of atoms 
unlike the clusters composed of 6-atoms. The RDM model of Phillips et al. leads to the 
existence of yN delocalized states in the linear chain lattice[12, 13]. In contrast, the ARDM 
model (for a^ = Eq) on a linear chain lattice leads to localized states only, with a localization 
length ^ given by [11] 

1.^ 4E^ _E^ p) 
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for weak disorder. Here E = 2 cos k is the energy band (in units of a constant nearest- 
neighbour hopping parameter V) of states of wavenumber k of the ordered lattice and .^o 
is the weak disorder Thouless localization length for uncorrelated Gaussian disorder in one 
dimension. The expression (2) is valid at energies E except for some special values within 
the energy band (including E = 0) where the weak disorder treatment must be corrected for 
Kappus-Wegner type anomalies[ll, 14]. It follows from (2) that the localization length near 
the band centre is strongly increased with respect to .^o by the effect of disorder correlation. 
For the tight-binding system on a ring threaded by a magnetic flux in the direction perpen- 
dicular to the ring the Schrodinger equation reduces to the set of difference equations[15, 16] 

- e ^ "t-O ipn+l - e ^ -fo ipn_i + En^n = Eipn, u = 2, 3, . . . , A^ - 1 , (3) 

•27r_^ - 271 ip 

— e ^ ■*() (y92 — e '^ ■*() (^f^ + Eii^i = Eifi , (4) 



e^'i'oipi-e '^ ■t-oipN_i + EN^N = EipN , (5) 



which embody the famihar flux- modified (twisted) boundary condition[3, 16] describing the 
effect of the magnetic fiux. Here 0o = hc/e is the flux quantum (with h the Planck constant, 
c the velocity of light and — e the electron charge), v^„ is the amplitude of an eigenstate 
wavefunction at site n, E and e„, are the corresponding eigenvalue and site energies in 
units of minus a constant nearest-neighbour hopping parameter. The equilibrium persistent 
current carried by the /c-th one-electron eigenstate of (3-5) is 

and the total persistent current in a ring containing < N^ < N spinless electrons in the 
zero temperature ground state is given by 

k 

where the summation extends over the A'^e lowest one-electron eigenstates of (3-5). 

Persistent current has been studied recently in the RDM of Dunlap et a/. [12] in several 
papers[17, 18, 19, 20]. In particular, Liu et al.[17] have shown that in the model of Dunlap 
et al. the persistent current is enhanced by the disorder correlation to values approaching 
the current in an ordered system if the Fermi level coincides with the energy of a delocalized 
state. On the other hand, their numerical calculations indicate that the persistent current 
is strongly suppressed when no delocalized state is present near the Fermi level [17]. 

The study of the persistent current in the ARDM is of interest in view of its speciflc local- 
ization properties discussed above. In particular this study may allow us to clarify whether 
the presence of a correlation induced delocalized state at the Fermi-level is indispensable 
for obtaining large persistent currents as discussed by the authors of Refs [17] and [18] for 
the RDM. These are our main motivations for the present work. The ensemble-averaged 
energy eigenvalues for the ARDM model for weak disorder are discussed in Sect. II, using 
the general perturbation expressions for the eigenvalues of Eqs. (3-5) derived in Ref. [16]. 
The detailed study of the ensemble-averaged persistent current successively for a system of 
Ne '^ N spinless electrons occupying a small portion of the energy-band of the ring de- 
scribed by (3-5), and for a system of N^. = N/2 electrons occupying the lower half of the 
energy-band at zero temperature, is presented and discussed in Sect. III. Some flnal remarks 
are given in Sect. IV. 



II. THE PERTURBED ENERGY LEVELS OF THE RING 

Exact general expressions for the one-electron energy levels of the ring threaded by a 
magnetic flux and perturbed by a weak site energy disorder have been derived from (3-5) 
in Ref. [16] , using an exact transfer matrix formalism and expanding the transfer matrix to 
second order in the disorder. In the absence of disorder the energy levels are given by the 
familiar tight-binding results[15, 16], 

El{<t.) = El = -2cos^-^(k + ^\ , (8) 

where /c = 0, ±1, ±2, . . ., and the solution of the eigenvalue equation yields the following 
exact expressions for the first- and second order perturbations of the energy levels [16] 

E, = El + EiHEfH... (9) 

^'^ ^ E''' = ^Y.^n , (10) 



n=l 



TV n-1 



Ef^ = ^ y TiE^'^ - Sm)iE^'^ - en) 

' iVsing,siniVg,A.A.^ -)^ -) ^^^^ 



X sin(n — 'm)gfcsin(A^ — n + m)qk 



where 




As recalled in Ref. [16] the weak disorder expansion breaks down at flux values equal to 
half-integer multiples of 0o including = where (11) diverges[21]. Such divergencies are 
familiar in other perturbative studies of persistent current [22] and are due to the degeneracies 
of the unperturbed electron energies (8) at these flux values. 

As in our earlier work[16] we are interested in the change of persistent current induced 
by the energy correction (11) averaged over the disorder. The derivation of the closed form 
of (11) averaged over the correlated disorder deflned by (1) involves the following steps: 

1. for each one of the four disorder factors (E^^^^, —€rnE^^\ E^^^Sn and €m£n (with E^^^ 
deflned by (10)) which enter in the site summations in (11), we identify explicitly 



the individual nearest-neighbour pair terms with a no n- vanishing correlation (1), in 
addition to the autocorrelation terms. 

2. after performing the disorder averaging, using (1), the various contributions reduce to 
combinations of geometric series over sites, which are then summed in closed form. As 
a guide to our calculations we list in the appendix the final explicit results obtained for 
the four distinct disorder terms in (11). The sum of the contributions (A2-A5) leads, 
after further reduction to the following exact expression for the averaged disorder effect 
in the ARDM: 



(4 = (^f)™^ + (4'r 
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(15) 



Here {El ''^'^n^orr rgprgggnts the averaged energy correction resulting from the random 
potential in the absence of correlation, which has been studied earlier[16], and {E^ 'y°'^^ 
represents the effect of the correlation of the disorder described by the dimer model of 
Ref. [11]. 

The equations (9) and (13-15) are used in the following section for studying the 
ensemble-averaged persistent currents. 



III. THE PERSISTENT CURRENT IN THE RING 

We first recall the results for the persistent current in the pure (non-disordered) tight- 
binding ring obtained from Eqs (6-8). The energy levels (8) occupied by the N^. spinless 
electrons differ for odd and for even N,,, as well as for different domains of relative magnetic 
flux 0/00) as is well known[15]. The sets of occupied /c-levels entering in the calculation 
of the total current in various cases may be inferred from the set of intersecting parabolas 
representing the energy level spectrum as a function of flux in the free electron limit (see e.g. 
flg. 2 in Ref. [15]). They are deflned by the values A; = 0, ±1, ±2, . . . , ±{Ne - l)/2 for odd 
iVe and -1/2 < 0/0o < 1/2; and by A; = 0, ±1, ... , ±{NJ2 - 1), -NJ2 for < 0/0o < 1/2 
and by A; = 0, ±1, ... , ±(iVe/2- 1), NJ2 for -1/2 < 0/0o < for even A^e, respectively. The 
summation of the contributions of the various levels (8) reduces to geometric progressions 
leading to the following mostly by well-known results for pure system currents /*^°-'[15, 17]: 



(0) sm(27r0/iV0o) 20 

-' ~ "-'0 '■ Tat — ~-'o- j~ 

sm -K /JS 00 

1 1 

odd iVe, -- < 0/00 < - 



(16) 



^(0) ^ ^^^ sin(7r/iV)(20/0o-l 
sin Tx /N 



■■ 777 - -^0 




even N,, < 0/0o < - ) , (17a) 

^(0) _ _j sin(7r/iV)(20/0o + l 




even iVe, -i < 0/00 < ) , (17b) 



J eVF 4:71 TlNe , . 

If' ^F = -^smA;j., kF = ^^ , (18) 

where the approximate limiting forms correspond to the limit N ^ oo and coincide in the 
case of (16) and (17a) with earlier analytical results for free electrons[15]. vp and kp in 
(18) denote the tight-binding Fermi velocity and the Fermi momentum, respectively[16]. 
The limiting free electron expressions in (16) and (17a, 17b) exhibit the familiar sawtoothed 
currents shown in flgs. (3. a) and (3.b) of Ref. [15] (whose legends must be interchanged!). 
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We now study the effect of the correlated disorder on the persistent current obtained 
from (6-10) and (13-15) for a system of N^ spinless electrons in the T = ground state. 
Our general discussion is made systematic by considering successively the case of low filling, 
Ne <^ A^, of the energy band and the case of higher filling exemplified by the half-filled band 
case, A^e = N/2. 

A. LOW FILLING OF ENERGY BAND 

For occupied levels with quantum numbers k <^ N and magnetic flux restricted by (16) 
and (17a, 17b) we use 

smq,c^q, = -^k + -j , (19) 

and by expanding (14-15) to lowest order for large N we obtain 

where 

(Ef)™ = iiLt^--^) , (21) 

^qk\ (po Nqk ' ' 
is the uncorrelated disorder perturbation at low band-filling obtained from (14). The Eq. 
(20) wil allow us to discuss the current generated by a system of A^e << -^ electrons. 

1. Single electron current 

It is instructive to briefly discuss the effect of correlated- and uncorrelated disorder on 
the persistent current for a single electron on the lowest level {k = 0) of the energy band (8) 
for -1/2 < 0/00 < 1/2. From (20) and (6) we obtain 



\-'fc=0/ disorder — -'0 | '7~ I ~T" 



where Iq = {47r^e)/Nh. 



^0 , el(f>o _27r(j) , , 0.2W Vo 



2n 



cot —— - 2(ei + 3a' 



sin2(27r0/0o) 27r0 0o ^° M 27r0 



, (22) 



It follows from (22) that the correlated disorder generally leads to an enhancement of persis- 
tent current with respect to the pure system result, Ik=o = — -^o20/A^0o, obtained from (8). 
In contrast the corresponding expression for the current induced by uncorrelated disorder 
vanishes for 0/0o ~^ and leads to suppression of persistent current at larger flux within 
the range (16). 

2. Low band filling 

At low band filling the occupied energy levels in the ring below the Fermi level corre- 
spond to small gfc- values (for large A^), so that the terms involving sing^ and related small 
quantities in denominators dominate in the disorder perturbation (14-15). By expanding 
the denominators in question to leading order in g^ and collecting terms we obtained the 
expression (20) for {E^, ) which we now use for calculating the persistent current induced by 
the system of iVg electrons in the ground state when 1 ^ A^e ^ A^- The averaged persistent 
current from the A;-th level is given by 



vre 1 

X-'fc/ disorder 



2hNqk 



e'oN , e^Q 27r0 2(e^ + 3a' 

COTj 



(23) 



sin'(27r0/0o) Qk 0o Nq^ 

A simple analytic procedure for performing the summation of the currents (23) from the A'e, 
lowest levels of the energy band (9) is to replace the summation by an integration. For this 
purpose we use the series summation formula of Euler-Maclaurin[23], 

T.fk= f ^ i\k)dk-\[i\K_) + f\K^)lK^ = k_-l,K+ = K + l , (24) 

where we have retained only the leading boundary term of the general formula[23]. Here A;_ 
and A;+ denote, respectively, the smallest and the largest value of the quantum number k 
labelling the occupied levels in the energy band of the ring. From the spectra of occupied 
energy levels detailed above for the domains —1/2 < 0/0o < 1/2 for odd and for even A'e 
we have: k+ = -k^ = {N^ - l)/2 for odd A'e; k+ = NJ2 - 1 and k_ = -Ne/2 for even A^e 
with > 0; and A;+ = A^e/2, k^ = —Ne/2 + 1 for even N^ and < 0. We thus obtain the 
following values for q^ at the lower and upper limits of the integral in (24): 
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(odd^..-l<^<i) . (25) 



«--^(-^= + f)- "--^(^= + ^ + 1) 



) 



(even iVe, -i < ^ < 0) (26a) 

2 00 

(even iVe, < ^ < ^) . (26b) 

00 2 

which are vahd for arbitrary < A^e < ^- Having discussed the case A'^g = 1 in III A 1 above 
we now focus mainly on values A'e ^ 1 for which 

gi.,^-gi._ = ^('l + 0(^i-) I, Ar,»i , (27) 

in the various cases detailed in (25-26). This will be sufficient indeed for illustrating the 
qualitative differences between these cases. 

The total persistent current (7) obtained from (23) involves three types of sums over 
occupied levels, namely 5*1 = J2k1k^^^2 = Y.klk'^ ^^^ ^3. = J2k1k^- ^^ recall 
the occupied levels in the various domains of flux, namely k = —{Ng — l)/2, — (A"e — 
l)/2 + 1, . . . , 0, 1, . . . (A^e - l)/2 for odd A^e and -1/2 < 0/0o < 1/2, k = -Ne/2 + 
1, -NJ2 + 2, . . . , 0, . . . NJ2 - 1, NJ2 for even A^^ and -1/2 < 0/0o < 0, and finally, 
k = -NJ2, -NJ2 + 1, . . . , 0, . . . NJ2 - 1, for even A^'e and < < 1/2. The series ^i 
which does not converge for k —^ oo requires special care. Using (12), we rewrite it in the 
following forms: 

^i = ^(l-2J), for odd AT,, -1<^<1 , (28) 

27r0 2 00 2 

'5i = ?^(l-2'/) + — , for even AT,, -i < A < q , (29a) 

27r0 Qk^ 2 00 

5i = ^(l-2J) + — , for even AT,, 0<^<i , (29b) 

27r0 qk_ 00 2 
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p I 



k=l 



(30) 



with p = /c+ in the case of (28), p = A;+ — 1 in the case of (29a) and p = k^ in the case of 
(29b). 

In applying (24) for performing the summation over the occupied levels in (30) we en- 
counter a difficulty which is the fact that the integral is not defined since its lower limit 
(zero) leads to an imaginary term. This difficulty may be resolved by introducing a small 
A;-cutoff for the domain of integration, say at k = kc = 1. This is justified since the lowest 
discrete levels contributing to the integral are the levels k = ±1, the effect of the level k = 
having been separated out in the first term on the right hand side of (28,29a,29b). With 
this regularization we obtain 



J 



In 



1- 



+ 



2(0g 



On the other hand by summing 5*2 and 5*3 using (24), we obtain 



«'iF) 



(31) 



N 



S, = — 



1 



2vr \ qk. 



qK. 



1 



nqt 



qL 



(32) 



S.^-^ 



Ik- 



(33) 



and using the definition (25) and (26a,26b) the leading explicit forms of (28), (29a,29b) and 
(32-33) are given by: 



Si 



iV0o 
27r(/) 



;i-2J), oddA^e 



(34) 



Si = ^il-2J)+ ^ 



Si 



2n<p 
27r0 



1-2J) 



-kN, 



N 



1 + 



1 + 



N, 



N, 



, even Nf. 



1 

- < — <0 , 

2 00 



, even N^-, < — < - 
00 2 



(35a) 
(35b) 



So 






^^N- 



(36) 
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s^^^li." 




(37) 



(po \7rNe 

where the hmiting forms (36) and (37) are vahd for even- as well as for odd Ne, in the 
domain -1/2 < 0/0o < 1/2. 

Finally, using (23) and (34-37), the change in the persistent current (7) due to the disorder 
is given by 



(/) 



loN' 



disorder 



^vrAT, 



2n 



L(i_2j)+g 



1 2TT(f) ^/ 1 



co.^.o(^) .m 



sin2(27r0/0o) T^^Ne 

where Q = for odd N^, Q = (ttN^)-^ for even N^ with -1/2 < 0/0o < and Q = -{ttN^)-^ 
for even A^e with < 0/0o < 1/2. 

It follows from (38) that for 1 ^ A'^e ^ ^, the effect of the disorder is dominated by 
the form of J in (31), which is negative for both signs of the magnetic flux 0. This shows 
that for odd A'^e the disorder reduces the persistent current with respect to the pure system 
result (16) while enhancing it with respect to the corresponding pure system values (17a) 
and (17b) for even Ng. 

We conclude that while the correlation effect may efficiently offset the effect of an uncor- 
related potential on the persistent current both in the case of a single electron (see IIIAl) 
and in that of a gas of a large number of electrons (see the case A^e = N/2 in IIIB below), it 
does not affect the effect of uncorrelated disorder on the persistent current in systems with 
1 << A^e << A^ at low order. 

B. HALF-FILLED BAND 

In the case of a half-filled band the terms proportional to 1/A^ and to 1/A^^ in the energy 
level perturbations (14-15) may be ignored relative to the remaining terms of order one for 
energy levels (8) near the Fermi level. In this case we thus have from (14-15) 

(4^y— ^-^cot^ , (39) 

'' 4smgfc 00 



(Ef V°- ^ y (cos g, - cot ^ sin q^) . (40) 
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The persistent current due to the effect of the uncorrelated random potential (39) has been 
discussed in Ref. [16] for a half-filled band for odd A'e in the range —1/2 < 0/0o < 1/2 
and for even N^. in the range < 0/0o < 1) using the Euler- Maclaurin series summation 
formula. However, throughout the present work we consider the typical domain —1/2 < 
0/00 < 1/2 for the detailed analyses of the effects of correlated disorder in the persistent 
current. This seems more useful since recent first numerical results for persistent currents 
in disordered rings with a correlated random potential have been presented for the range 
— 1/2 < 0/00 < 1/2 only[17, 18, 19]. From (39) the effect of the uncorrelated random 
potential on the persistent current in the /c-th level is given by 

(4)unco.. = :^ ■ ,'\,,. ■ +Q(l/iV) • (41) 

Ih sm (27r0/0o) sm q^ 

whose summation for the N^ lowest occupied levels of the energy band (8) in the domains 
— 1/2 < 0/00 < and < 0/0o < 1/2, respectively is performed using the Euler- Maclaurin 
formula (24). Using (18), this yields the following final result for the effect of the uncorrelated 
disorder on the mean persistent current in the half-filled band: 



(h - I ^'° 

\ -* / uncorr -* , 



;sin^(27r0/0o) 



A^ (l + cosgx+)(l -cosg/^_; 

In ; 

A-K [1 — cos qKj^)[^ + cos qK_, 




, (42) 

in terms of the g^- values (25,26a,26b) for iVe = N/2, near the upper and lower boundary 
values A;+ and A;_ of the domain of occupied A;-levels and for flux values within the interval 
(— 0o/2, 0o/2). The equation (42) may be readily simplified for Ne = N/2 ^ 1 by expanding 
the square bracket to leading order in 1/A^, using (25, 26a, 26b). This yields 

8sm (27r0/0o) 
A = ^ + 0{1/N'), oddiV,, -1:<^<1 , (44) 
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A = 1 + ^ + 0(1/Ar2), g^g^ ^^^ _1 < A < 0, , (45a) 

00 2 00 

A = -l + ^ + 0(l/iV2), eveniV^, < ^ < ^, , (45b) 

00 00 2 

where (43) with A defined by (44) or (45b) has been obtained earher in Ref. [16]. 

The persistent current due to the correlation for an electron on the k-th level is, from (6) 
and (40), 

2 

716 a 

(/fc)corr = - — — iTT^TTTTTsingfc + 0{l/N) , (46) 



h sin (27r 

which leads to the following exact results for the persistent current for the system of N^ 
spinless electrons occupying the lower half of the energy band of the system: 



IT\ J ^^' sin(27r0/Ar0o) 

\-' /corr ^ _. 2/ 



4sin^(27r0/0o) simi/N 



-/« m ""' 



\-' /corr -'( 



V0O/ 4sin^(27r(; 
(odd iV„ -1/2 < 0/00 < 1/2) , (47) 



Na^ sinffgil 



°4sin2(27r0/0o) sirni/N 
No" 



^°Uo ^V 4sin2(27r0/0o) 
(even iVe, +sign for - 1/2 < 0/0o < 0, -sign for < 0/0o < 1/2) , (48) 

By comparing these expressions with the pure system currents (16) and (17a, 17b) it follows 
that the effect of the correlation of the random potential is, in all cases, to enhance the 
persistent current. 

Finally the total average persistent current (/) in the ring is given by the sum of free 
electron currents [large A^-limits of pure tight-binding system currents, Eqs. (16) (17a) (17b)] 
and the changes in the persistent current induced by the disorder via the correlations of the 
random potential [Eqs. (47-48)] and via the random potential in the absence of correlation 
[Eqs. (43-45)]: 
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(/) 



Njel - 2a^) 
^sin2(27r0/0o) 



11 Ar 1 1 

.odd A^e, — < — < - 
2 00 2 



(49) 






-^°'^^^ 
-I- 



— ^^ -^ , even A^^, — < — < 

^sin^ 27r0/0o y 2 0o 



Njel - 2q^) 
8sin2(27r0/0o^ 



1 
,even iVp,0 < -- < - 



(50a) 
(50b) 



with Iq defined in (18). These expressions are invahd near = 0, which is a consequence 
of the breakdown of the perturbation expression (11) for flux values equal to half- integer 
multiples of the flux quantum, as mentioned earlier. 

The expressions (49) and (50a-50b) indicate that the disorder acts to renormalize the free 
particle persistent currents (16) and (17a-17b by a common factor. In particular, it follows 
that the correlation of the random potential leads to a systematic increase of the persistent 
current proportional to a"^. However, the important feature of (50a-50b) is that for a"^ less 



than the critical value a^ 



a^ = el/2 the persistent current is reduced by the correlated 
disorder while for a^ > a"^ it is enhanced by it for any parity. In this respect, the results for 
the persistent current at low band fllling, 1 <^ N^, 'ti N, discussed in 111 A above show an 
interesting asymmetry. Indeed, in this case the effect of the correlated disorder is to reduce 
the persistent current for odd A'e and to enhance it for even A'e. This is a consequence of 
the fact that the flux dependence of the persistent current (38) induced by the disorder does 
not follow the simple linear scalings of the pure system currents (16, 17a, 17b). 

C. DISCUSSION OF THE RESULTS 

In this subsection we discuss the restrictions on ring lengths A^ and on relative flux values 
in (49) a nd (50a-50b) which are imposed by our non-degenerate perturbation study of the 
effect of the disorder on the energy levels of the ring in the presence of the applied flux. 
We focus mainly on the case of the half-fllled energy band where the effects of the disorder 
on the persistent current may be observed in wide range of nanoscale ring parameters, as 
shown below. 

First, the perturbation expansion (9) in terms of the random site energies relative to the 
flxed intersite hopping parameter implies, by deflnition. 
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|£o| << 1, |a| << 1 , (51) 

which reflects the fact that the first order correction to the energy levels, E^. , of the pure 
system is small. Next, the typical value of the first order shift of energy levels due to the 
disorder, \l {Ef, ) = {e1 + a^)/N must be small compared to the spacing, |A|, of the energy 
levels of the pure system. From (8) we obtain 







which leads to the condition 



VnJeI + a2 << 47r sin -^ ( - + y- ) ~ 47r , 



(52) 



at the Fermi level (where k = ±Ne/2) in the half-filled band case. This condition which is 
important for close-spaced level systems (large N) expresses that disorder acts as a small 
perturbation on the (observable) excitation frequencies defined by A. Since, for band fillings 
A'^e << A^, the Eq. (52) restricts the discussion to rings of perimeters less than a few 
nanometers, we focus, in the remainder of this subsection, on the half-filled band case 
iVg = N/2, where our analysis covers a much wider range of nano metric rings sizes. Also since 
here we are mainly interested in analyzing the effect of small random potential correlations 
(a^ < £g) for rings of fixed length A^ we replace the condition (52) by the slightly weaker 
form: 



^N^Jel « 4:71 (52.a) 

at the fermi level of the half-filled band. 

Finally, an important restriction on the allowed flux values in (49) and (50a), (50b) arises 
from the second order perturbation expression for the effect of the disorder in the energy 
levels (9-11). This requires the typical second order corrections of the energy levels to have 
magnitudes of the order of Eq. In the half-filled band case this condition is generally obeyed 
by Eqs (39-40) for {E^, ) at the fermi level (where sing^ ^ 1) by requiring cot -r^ ~ 1. This 
yields a lower limit 

0mm = 0, 12500 , (53) 
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for the flux values at which the perturbation treatment of the disorder is valid. 

We close this subsection with some typical numerical results for the renormalization factor 
in (49), 50a and (50b), 

^"^ Ssin^^ '"^-el ' ^^^^ 

<P0 " 

of the pure system current by the effect of the disorder as a function of the correlation param- 
eter 1] = a^ /el. We choose |eo| = 0.1 and define a typical ring length A^ = (0,47r)^/eQ = 158 
compatible with (52. a). Using these values, we obtain the results listed in table 1, for var- 
ious values of r], for = ipmin and for an intermediate value, = 0.2500) in the domain 
< 0/00 < 1/2. Note that increasing the flux in the case r] < 0.5 enhances the magnitude 
of the disorder effect on the persistent current, while in the case rj > 0.5, increasing the flux 
reduces the effect of the correlated disorder. 

IV. CONCLUDING REMARKS 

The main conclusions of this work have been discussed in Sect. Ill along with our detailed 
results for persistent currents in a metal ring in the presence of correlated disorder and the 
discussion of their domain of validity. 

We conclude with some further general remarks, particularly in connection with the yet 
unexplained experimental observations of persistent currents. 

By studying the persistent current in sufficiently large ensembles of rings it has been 
possible recently to obtain the correct sign, in addition to the magnitude, of the current [24, 
25, 26]. These measurements have shown the persistent current to be diamagnetic at low 
fields (i.e. -^ < 0, ^ 0). This observation constitutes a severe additional constraint for 
the validity of theoretical explanations of persistent current in disordered rings. Indeed, as 
shown by (16), (17a) and (17b), a diamagnetic persistent current for any parity of the number 
of conduction electrons is generally obtained only for non-interacting free electrons. Now, our 
analytical expressions (49), (50a) and (50b) for the effect of disorder on the persistent current 
in a half-filled band includes an additional flux- dependence which prevents the current from 
having a flxed negative sign independent of parity and ring size. We observe, however, that 
our theoretical results for persistent currents are inapplicable to measurements in rings of 
mesoscale perimeters of 2 to 8 /xm used in the various experimental studies[4, 5, 6] and 
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[24, 25, 26]. Indeed, for values I^qI — 0.1, whith corresponding values < |q;| < 0.1, the 
condition (52) restricts the validity of our results to rings with a few hundred atomic sites 
i.e. to nanoscale perimeters less than 100 nm. We note incidently that typical ring sizes for 
which the effect of correlations of the random potential has been analyzed numerically in 
Refs.[17, 18, 19] lie also in this range. In conclusion, it would be interesting if experimental 
studies could be performed for ensembles of rings of nanoscale sizes in order to test the 
detailed behavior of persistent currents shown theoretically in this paper and numerically, 
for a different model, in Refs.[17, 18, 19]. 

Our results for the persistent current for a half-filled band in our ARDM model indicate 
a strong enhancement of the current due to the correlation of the random potential for 
o? > ^o/2- The detailed mechanism for the enhancement of persistent current by correlation 
for the ARDM model [11] on a ring threaded by a magnetic flux is quite different from the 
mechanism identified in recent numerical calculations for the RDM model of Dunlap, Wu and 
Phillips[12] in Refs. [17, 19]. We recall that the linear chain RDM model differs essentially 
from the ARDM by the existence of a fraction \fN of extended states[12, 13]. For the 
RDM model values of persistent current of the order of the values for free electrons are 
found only when the Fermi level coincides with the energy of an extended state, both for 
even A^e[17, 19] and for odd A'^e (for an asymmetric dimer model) [18]. In contrast, for N^ 
values corresponding to relatively large band fillings but for which the Fermi level does not 
coincide with an extended state, the persistent current is found to be strongly diminished 
by the correlated disorder [17, 18]. Thus for the RDM a strong compensation of the current 
reduction due to the uncorrelated disorder by the effect of correlation does not seem to 
exist [17, 18], in contrast to what is shown for the half-filled band in (49) and (50a,50b) for 
the ARDM. We note, however, that the numerical results of Refs[17, 18] pertain to a case 
of maximal disorder (equal number of a- and 6-atoms in the RDM), while our analytical 
results are restricted to weak disorder in the ARDM. 

As shown in Sect. Ill C, our results for the persistent current are valid only for flux values 
larger than a value of the order of 0o/8- Thus the sign of the persistent current obtained form 
Eqs. (49, 50a, 50b) could be compared with experimental results at flux values > 0o/8 
only, rather than for ^ as in recent experimental studies[24, 25, 26]. 
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V. APPENDIX 



Defining 



/i„ 



sin(n — m)gfc sin(A^ — n + m)qk 



sin Qk sin A^ qk 

(2) 

the various contributions to the disorder average, {El ), of (11) are given by 



^ N n-1 



(eg + a') 



N 



{N cot Nqk - cot gfc) , 



n=2 m=l 



-. Af n-1 



N 



n=2 m=l 



AN sin qk 

a^ sm{N — l)gfc 
2 sm qk 



(A.l) 



(A.2) 



(A.3) 



^ Af n-1 
n=2 m=l 

-2 



(^g - «^) 



(A^cot A^gfe-cot qk) + 



a 



4A^singjfc^ ^"' ^""^ ' 4A^^ sin g^ sin A^g^ 

g-j(Af+2)gfe /e*^*^'=(l — e^*^*'=) A^\ 
1 _ e-4«9fc ( iZT^k y j 



N^ cos A^gfc 



-2 



1 _ e-4i<?fe \ 1 — e^*'"= 



c.c. 



, (A.4) 



N n-1 



i5:5:{B'"oA 



(A^ cot Nqk-cot qk) + 



n=2 ni=l 



4A^2 sin gfc sin A^g^ 1 

-i{N-2)qk /]_ _ g2i(Af-l)(3fc ]^ _|_ ^-2iqk _ ^2iNqkl^Aiqk _|_ g-2j(7fe\ 



1^5 - « J 
4A^ sin gfc 



2A^(A^-l)cosA^gfcl 



1 _ Q-^iqk \ X _ e2ji3fc 



- (A^ - l)e'*^'= - — + 1 e 



A^ 



X _ g4i<7fe 



= -2iijfc Q-6iqk 



C.C. 



. (A.5) 



The sum of (A.2-A.5) leads after appropriate reductions to the final form in (13-15). 
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V 


4>o 


A 





0,125 
0,25 


0,605 
0,803 


0,2 


0,125 
0,25 


0,769 
0,882 


0,4 


0.125 
0.25 


0,921 
0,961 


0,5 


- 


1 


0,6 


0,125 
0,25 


1,079 
1,040 



TABLE L Disorder dependent free electron persistent current renormalization factor (53) for var- 
ious values of the disorder correlation parameter r/ 
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